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ABSTRACT 

It is argued that the superfluid core of a neutron star super-rotates relative to the 
crust, because stratification prevents the core from responding to the electromag- 
netic braking torque, until the relevant dissipative (viscous or Eddington- Sweet) 
time-scale, which can exceed ~ 10 3 yr and is much longer than the Ekman time- 
scale, has elapsed. Hence, in some young pulsars, the rotation of the core today is 
a fossil record of its rotation at birth, provided that magnetic crust-core coupling 
is inhibited, e.g. by buoyancy, field-line topology, or the presence of uncondensed 
neutral components in the superfluid. Persistent core super-rotation alters our 
picture of neutron stars in several ways, allowing for magnetic field generation 
by ongoing dynamo action and enhanced gravitational wave emission from hy- 
drodynamic instabilites. 

Subject headings: dense matter — hydrodynamics — gravitational waves — stars: 
interiors — stars: neutron — stars: rotation 



Introduction 



Shear flows are endemic to neutron stars. For example, the interpenetrating com- 
ponents of the stellar interior rotate differentially. In the inner crust, an inviscid neu- 
tron condens ate flows through a rigid nuclear lattice. Nucle ar pinning controls the rela- 
tive velocity (IHaskell et al.ll2011r IWarszawski fc Melatosll201ll ). The condensate also flows 
relative to a viscous fluid component, consisting of electrons, protons, and uncondensed 
neutrons. Entrainment, hydromagnetic stresses, and Ekman pumping control the rela- 
tive velocity; under certain co nditions, the nucl e ar lattice and vi s cous component lock 



together magnetically in < 1 s (jAlpar et al.l Il984 : Ijahan-Miril 120021 : ISidery fc Alparl 12009 
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van Eysden fc Melatoal2010l ). An gular momentum store d in the condensate is released errati- 
cally to the cru st, causing glitches (phemar fc Lvndll996f). although the trigger mechanism re- 



mains unclear ( Warszawsk i fc Melatos 



Warszawski fc Melatos! 120 111 ) 



2008: 



Melatos fc Warszawskill2009l ; ISiderv et al.ll2010 



Additionally, differential rotation can occur between spatially separate regions inside the 
star. Examples include inductive and capacitive zone s, invoked to exp l ain the wide rang e 



of post-glitch recovery time-scales in the same pulsar (lAlpar et al. 



core is a crystalline color superconductor ( lAlford et al. 



1996 



2007l : lHaskell et al-koill ): and the crust a nd core JPeralta et al.l Eo05l.l2006bh. especially if the 



Wong et all 12 001); 



regions of strong and wea k pinning, predicted by nuclear theory ( Jones! 1991 ; Avogadro et al. 



20081 ) . Other condensed astronomical 



bodies are known or suspected to sustain crust-core shear, notably the Earth, where core 
super-rotation has been observed in seismi c travel-time data, and the size of the effec t and 
its variability time-scale are hotl y debated (pong fc Richards! 1 19961 ; I Waszek et al.ll201ll ). and 
the planets Jupiter and Saturn (jGuillotlll999l ). 



One natural cause of crust-core shear is gravitational stratification. Oceanographic and 
other studies of Ekman pumping show that stratification limits the volume of fluid that spins 



up (down) in an impulsively accelerated (decelerated) container (|Walin 



Ekman time-scale. 



St-Maurice &; Veronisl Il975l ; lAbney &; Epstein! Il996l ; iDuck Foster! l200ll ) , at least over the 



1969; Clark 1973 



The entire fluid does corotate eventually on a much longer, dissipative 



time-scale ( [Greenspan fc Howard! Il963l ; IClarklll973l ). after inertial oscillations die away. But 
by then, in a neutron star, the crust has decelerated further electromagnetically, so a crust- 
core shear always pers ists. The shear shapes the gravitational wave sig nal emitted during 
post-glitch relaxation (Ivan Eysden fc Melatos! 120081 ; Bennett et al.ll2010l ). 



In this paper, we apply classic ideas from Solar spin down ( jSakurai et al.lll971t IClark 



19731 ) to show that the core of a neutron star rotates systematically faster than the crust in 
response to the combined action of stratification-limited Ekman pumping and electromag- 
netic braking. In §|2j we identify the processes that govern the approach to a steady state in 
stratified spin down. In ^3], we evaluate the relevant time-scales and estimate the crust-core 
lag, finding that it is large early in the star's life but decreases with age. In §U we ask 
whether magnetic crust-core coupling is strong enough to erase the crust-core lag and, if so, 
under what conditions, a subtle question. The observational consequences of a persistent 
lag for precession, magnetic field evolution, and gravitational radiation are outlined briefly 
in £j5j Each consequence deserves to be assessed in detail, a challenging task which awaits 



^n this paper, we restrict attention to the spin-down problem relevant to an electromagnetically braking 
neutron star. 
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further work. 



2. Stratified Ekman pumping 



The fluid interior of a neutron star responds to changes in the angular velocity Q of the 
crust via Ekman pumping. In an unstratified, Newtonian fluid, Ekman pumping operates 
as follows: the Coriolis force drives a secondary circulation in boundary layers of thickness 
E X I 2 R abutting the crust; then the boundary flow sucks flui d from the interior, spins i t 
down, and recycles it on the Ekman time-scale = E~ l l 2 VL~ l ( Greenspan fc Howardlll963l ). 
The Ekman number, E = u/(R 2 Q), is a function of the kinematic viscosity v and stellar 
radius R. It is the only dimensionless variable in the problem. In an unstratified, two- 
component, Hall-Vinen-Bekharevich-Khalatnikov superfluid, much the same thing happens 
to the v iscous component, which in turn drags along the inviscid component via mutual 
friction (Ivan Evsden fc Melatoslboiol . boill . l2012h . 



Stratification alters Ekman pumping fundamentally. In addition to E, two new dimen- 
sionless variables enter the problem: the Burger number, S = N 2 /Q 2 , which quantifies the 
relative strength of the buoyancy and Coriolis forces, and the Prandtl number, a = v/k, 
which quantifies the relative rate at which momentum and heat are transported diffusively. 
The symbols N and k stand for the Brunt- Vaisala frequ ency and the r mal d iffusivity respec- 
tively. We see below that, in a neutron star, S is large [cf. iBaym et al.l (119761 )]. and a is either 
small or of order unity, depending on exactly how superflui dity affects k. The interplay and 



ordering of E, S, and a lead to a rich variety of motions (ISakurai et al.lll97ll ; IClarkl Il973 
St-Maurice fc Veronislll975l ; ISpence et al.lll992l ). 



First and foremost, stratification limits the interior volume influence d by Ekman pump- 
ing t o a boundary region, called a buoyancy layer, of thickness S~ l l 2 R ( Walinlll969l ; I Clark 
19731 ) . In a neutron star, the buoyancy layer is thicker than the unstratified Ekman layer 
(S <C E^ 1 ), but it does not suck in, spin down, and recycle interior fluid to produce uni- 
form rotation. Instead, on the Ekman-like time-scale S~ 1 / 2 E~ 1 / 2 £l~ 1 , the interior organizes 
itself into a quasi-steady, sheared, azimuthal flow, whose angular velocity runs smoothly 
from its original value deep in the core t o match the spun-down buoyancy layer at the crust 
(ISakurai et al.lll97ll ; ISpence et al.lll992l ). The angular velocity profile is set by the geome- 
try an d temperature boundary conditions; in general , Taylor-Proudman columnarity breaks 
down (jSt-Maurice &: Veronislll975l ; lFriedlanderlll976l ). 



Secondly, stratification distorts the boundary layer (lWalinlll969l ). In a cylinder, the side 
wall cannot carry a vertical (i.e. falling) mass flux, so it cannot accept an 0(E 1 ^ 2 ) Ekman 
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inflow. As a result, fluid is pushed in to the comers and thence erupts into the interior as a 
jet, leaving behind a stagnant region (jSpence et al.lll992l ). Er uption also occurs in a sphere, 
with Ekman suction turning into blowing at latitudes < 24° ( jClark et al.lll971t iFriedlander 
19761 ). The jet exists for Scr > E 2 / 3 , which always holds in a neutron star. 



Ultimately, of course, the sheared interior flow must evolve towards uniform rotation. 
When does this happen? Not until the faster of two dissipative processes takes hold: viscous 
diffusion (viscous time-scale t u = E^ 1 ^ 1 ), or thermally driven circulation (Eddi ngton-Sweet 
time-scale t K = aSE~ l Vt~ l , i .e. the thermal di ffusion time-scale multiplied by S) ( 1Clarklll973l ; 



St- Maurice &: Veronisl 1 19751 ; iHollerbachl 12009 ) . If thermal conduction in the neutron super- 
fluid is mainly collisional, one has a ~ 1 (IBaiko et al.l l2001t lUrpinl 120031 ). and it is safe 
to ignore the latter process. However, in terrestrial superfluids like helium II, convective 
counterflow of the small fraction of uncondensed neutrons boosts the effective thermal con- 



ductivity dramatically (ICutler fc Lindbloml 119871 ; IBaiko et al.l 120011 ) . The possibility that 
similar physics is at work in a neutron star is hard to quantify but cannot be dismissed, 
so we analyse it too in what follows. For aS < 1, a thermal layer diffuses inwards, its 
thickness growing with time t as (nt) 1//2 . It encounters the inner edge of the buoyancy layer 
at t — crS' _1 -E ,_1 fi, whereupon the angular velocity decrement imposed by the crust diffuses 
i nwards within a layer whose thickness grows oc t 1 / 4 , reachi ng the center once t equals t K 
JSakurai et al.lll97ll : Iciarklll973l : Istldaurice fc Veronisl 1 19 75^ . 



3. Long-term crust-core shear 



The maximum angular velocity gradient within the star depends on where the electro- 
magnetic braking time-scale Q/Cl fits into the hierarchy of hydrodynamic time-scales. Over 
the long term, there are two astrophysically interesting scenarios: £g < Q/Cl < min(t u ,t K ) 
and mm(t u ,t K ) < Q/Q. In the first scenario, early in the star's life, the angular velocity gra- 
dient equals (f2 — &)/R (if t v < t K ) or (fl — ^)(^/^k) 1//4 /-^ ^ > ^k)j where t is the star's 
age, and Qq is the initial angular velocity; in other words, the rotation of the core is a fossil 
record of its rotation at birth. In the second scenario, when the star is older, the angular 
velocity gradient equals |0| mm(t u ,t K )/R; all memory of the natal rotation disappears, and 
the gradient is proportional to the s mall amount of spin d own that o ccurs while t he relevant 
dissipative process unfolds [see §6 in Sakurai et al. ( 1971 ) and §5 in Clark ( 19731 )]. § 



2 For aS -C 1. iFriedlanden (|1976l ) found that uniform rotation is es tablished on the visc ous time-scale t v , 
if the crust experiences a constant spin-down torque, contrasting with lSakurai et al.l (|197ll ). who found that 



t K is the corotation time-scale, if the crust decelerates impulsively. This subtlety is an idiosyncracy of the 
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To quantify the shear, we refer to the latest, ab initio calculations of transport coeffi- 
cient s in bulk nuclear matter ( lYakovlev et al.l Il999l ; iBaiko et al.l l200ll ; IShternin fc Yakovlev 
20081 ) to evaluate the dimensionless hydrodynamic variables for a standard neutron star with 
mass w 1.4M and R = 10 6 cm: 
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In equations (P)-©, p and T denote the density and temperature in the core, p — 2.8 x 
10 14 gcm -3 is the nuclear saturation density, and R v and R K are superfluidity suppression 



factors for v and k respectively. □ □ We assume npe matter for simplicity: the electrons are 
i deal, degenerate, and ultrarelativistic, while the nucleons are nonideal and nonrelativistic 
( lYakovlev et al.l Il999l ). To calculate the thermal conductiv ity, we assume th at neutron- 
neutron and neutron-proton quasiparticle collisions dominate ( IBaiko et al.ll200lf ) (but set the 
dressed masses equal to the bare masses for simplicity), neglect electron conduction, neglect 
superfluidity in the scattering frequencies (but not elsewhere), and calibrate the result against 
Figure 1 in the latter reference. To calculate the shear viscosity, we assume that electron- 
electron collisions dominate, except at T ~ 10 7 K where neutrons also contribute, we neglect 
muons (consistent with the npe assumption), and we include th e important effect of Landa u 
damping mediated by weakly screened transverse plasmons ( IShternin fc Yakovlev! 120081 ). 
Calibrating against Figure 4 in the latter reference, we see that Landa u damping lowers v 
by ne arly an order of magnitude with respect to the traditional value ( jCutler fc Lindblom 
19871 ). and v scales as T -5 / 3 instead of T -2 , which is astrophysically significant because it 



regime aS <C 1, which is unlikely to be reached in a neutron star, even if superfluidity enhances k. 

3 In general, R v and R K are complicated functions of T, with R v ^ R K . As a rough rule of thumb, one 
finds R v = R K = 1 in the absence of superfluidity, R u « R K ps 0.2 if only the neutrons are strongly superfiuid, 
and R v ~ R K ~ 0.05 if both the neutrons and protons are strongly superfiuid. If convective counterflow in 



the uncondensed neutrons facilitates thermal transport, the foregoin g approximations still apply to R, 



we have R K ^> R v , but the value of R K jR v is not known accurately (jYakovlev et al 



1999: 



Baiko etal 



and 



2001 



Shternin fc Yakovlevll2008l) . 



4 The core temperatu re is related to the surface temperature via a two-zone heat-blanket model 
(jGudmundsson et al.lll982l ). and the surface temperature is related to the age t via 'fast' (superfiuid) cooling 
curves (including kaon and pion condensates and the di rect U r ea pro cess) calibrated against 11 neutron stars 
with measured thermal X-ray emission [see Figure 2 in IPage (1995 
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leaves a T dependence in a. The Brunt- Vaisala frequency, arising from departures from beta 
equil ibrium, is proportional to the s quare root of the electron-neutron ratio and is taken 



from iReisenegger & Goldreichl (119921 ). 



Figure [T] illustrates how much core super-rotation to expect in neutron stars of different 
ages. The results are expressed in terms of the Rossby number, Ro = AQ/Q, where AQ 
is the angular velocity of the core relative to the crust. One obtains Ro = Q /Q — 1 and 
Ro = |0| mm(t u , t K )/fl respectively in the two scenarios described at the start of §3J Contours 
of Ro are graphed on the plane spanned by the rotation period P and its time derivative P. 
Known radio pulsars are overplotted as blue dots. The blue band at the left of the diagram 
marks the excluded region S < 1, where stratification is weak, buoyancy does not impede 
Ekman pumping, there is no fossil rotation, the long-term shear Ro ~ tg/t < 10 -6 is too 
small to be of interest, and the analysis in this paper does not apply. 

Consider first the example R v = R K = 0.05 (red, dashed contours), which corresponds 
to strong superfluidity in the neutrons and protons, without counterflow-enhanced thermal 
conduction. In this regime, we find t v < t K ; viscous diffusion pre-empts thermally driven 
circulation throughout the displayed subset of the P-P plane, and the orange shading in 
the figure (explained in the next paragraph) is irrelevant. In most of the figure, below and 
to the right of the green shading, we also have t > t v and hence Ro = t u /(2t). The red 
contours run from Ro = 10 _1 to 10~ 6 , spaced by one dex. Numerous objects with spin-down 
ages < 10 4 yr are predicted to have relatively high Rossby numbers. In particular, there 
are three objects with Ro > 10 _1 (the young pulsar PSR J1846— 0258 in the supernova 
remnant Kesteven 75, and the magnetars PSR J1550— 5418 and PSR J1907+0919), and 18 
with 10~ 2 < Ro < 10" 1 (including three with Ro « 0.1). PSR J1846-0258, which is only 
0.7 kyr old, and whose braking index has be en measured by abs o lute pulse numbering;, i s 



an especially good rotation fossil candidate (IGotthelf et al.l l2000t iLivingstone et al.ll2006l ). 



In the top left corner of the figure, outlined by the green triangle, where young pulsars 
with t < 10 3 yr and T > 10 8 K reside, we find t <t v and hence predict Ro > 1, i.e. the core 
retains much of its natal angular velocity. No radio pulsars have been discovered in the green 
triangle to date, but such objects, whether young like PSR J1846— 0258 or magnet ar- like, 
are easily accommodated within our current picture of neutron star evolution. 

Next consider the example R v = 0.05, R K = 5 x 10 2 (black, unbroken contours), which 
models roughly the counterflow effect. The orange shading now becomes relevant; it divides 
the P-P plane into two parts. In the white region to the right of the orange curve, we 
find t v < t K , and nothing changes from the previous example: viscous diffusion dominates, 
and the black Ro contours overlap the red ones. To the left of the orange curve, however, 
the opposite is true: we find t K < t u , and thermally induced circulation occurs faster than 
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Fig. 1. — Contour plot of Rossby number as a function of spin period P (horizontal log 
axis, in s) and period derivative P (vertical log axis, in ss _1 ), running from Ro = 0.1 (top 
contour) to 10~ 6 (bottom contour) in steps of 1 dex. Red, dashed curves: strong neutron 
and proton superfluidity, collisional thermal conduction (R u = R K = 0.05). Black, unbroken 
curves: counterflow-enhanced thermal conduction (R u = 0.05, R K = 5 x 10 2 ). Blue dots 
mark known radio pulsars. The colored regions cover the following regimes. Blue: S < 20, 
stratification too weak to impede Ekman pumping, Ro negligible. Heavy green: fossil stage, 
t < t u < t K . Light green: fossil stage, t < t v < t K (t < t K < t u ) for R K = 0.05 (5 x 10 2 ). 
Orange: t > min(t u ,t K ) = t K , natal rotation erased, Ro decreases monotonically with age, 
thermally driven circulation pre-empts viscous diffusion. White: t > min(t^, t K ) = t u , viscous 
diffusion pre-emp ts the r mally driven circulation. Other parameters: p = po, Qq = 200n s _1 
(fiducial), T from IPagd (1l998[ ) (see footnote BJ. 
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viscous diffusion. Below the irregularly shaped, heavy green region, we also have t > t K — 
9.7(RJ5 x 10 2 )- 1 (p/p ) -1/3 (T/10 8 K) 2 (n/10 2 s- 1 )- 2 yr and hence Ro = t K /(2t). That is, 
when R K is large, the green region to the right of the orange curve (t u < t K ) is a promising 
hunting ground for rotation fossils with Ro > 1; we do not expect to find many to the left of 
the orange curve. Likewise, when R K is large, it is harder to find objects with Ro less than 
unity but still relatively large; as the black contours tilt upwards in the orange region, one 
needs a younger, faster pulsar to achieve a given Ro than when R K is small, yet young pulsars 
spin down quickly. A meaningful fraction of the pulsar population lies within the orange 
region. No known objects are predicted to have Ro > 0.1, three have 1CT 2 < Ro < 1CP 1 [the 
gamma-ray pulsar PS R J1023— 5746, discovered in a blind Fermi Large Area Telescope search 



flSaz Parkinson! l2010f ). PSR J0540-6919 in the Large Magellanic Cloud, and the Crab], and 
16 have 10" 3 < Ro < 1(T 2 

Figure [1] also maps the history of differential rotation in individual objects. As a pulsar 
brakes electromagnetically, it moves from the top left to the bottom right of the figure along a 
diagonal spin-down trajectory, PP = constant. If the trajectory intersects the green triangle, 
the early part of the star's life is spent as a rotation fossil. □ When the star crosses from 
the green region into the orange or unshaded regions, Ro drops suddenly on the time-scale 
mm(t u ,t K ), as dissipative processes take over. Subsequently, the crust-core shear decreases 
monotonically as the star brakes and cools, with Q/Q oc t^ 1 and hence Ro oc t _1 T 5//3 or 
Ro oc T 2 in the unshaded or orange regions respectively. § If the trajectory starts outside 
the green triangle, the latter behavior characterizes the star's whole life: there is no fossil 
stage. A wide range of terminal Ro values are po ssible, when the trajectory intersects the 



pair-cascade death line (IHibschman fc Aronsll200ll ). 



4. Magnetic coupling 



How does the stellar magnetic field alter the spin-down scenario above? This question 
turns out to be subtle. The physics depends sensitively on the internal magnetic topology 
and the multi-component nature of the superfluid interior, both of which remain highly 
uncertain. A proper understanding also requires a theory of Ekman pumping in a stratified, 
magnetized fluid in a sphere, which does not exist yet in the literature even for uniform 
magnetization, let alone for a realistic (e.g. linked poloidal-toroidal) field geometry or for a 



5 For the minority of neutron stars born in the blue region, where stratification is weak, Ekman pumping 
erases the natal rotation on the fast time-scale ts- 

6 The transition between the two scalings occurs smoothly at the wavy orange curve in Figure [l] 
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superconductor. In this section, therefore, we restrict ourselves to canvassing some of the 
main issues and presenting associated order-of-magnitude scalings, noting that large-scale 
numerical simulations will ultimately be required to settle the matter decisively. 

Magnetic coupling in its simplest guise (i.e. in an unstratified, uniformly magnetized 
star) speeds up angular momentum transport dramatically, enforcing corotation between the 
crust and core over ~ 10 2 s for typical magnetic field strengths and spins. The process occurs 
in two stages. Firstly, the charged component of the superfluid is accelerated by Lorentz 
forces, mediated not by Alfven waves but rather by rotation-dom inated, hydromagnetic- 
inertial waves, if the star is not superconducting ( lEassonl Il979bf ). or by cyclotron- vortex 
waves, whose r estoring force comes from the vortex tension, if the star is superconducting 
( IMendelll Il998[ ). Secondly, the charged fluid drags along the neutron condensate via mutual 
friction, as cha rges scatter off ne utron vortices magnetized by the Fermi liquid interaction 
(entrainment) (lAlpar et al.lll984j ). Overall, for uniform magnetization, the coupling time- 
scale is given by t B = 2(V A /Ml)- 2 ^- 1 = 6.6(p/p ) 1/3 (^/10 2 s~ 1 )~ 1/3 (5/10 12 G)~ 2/3 s in 
a non-superconducting star ( lEassonl 1 19 79b| ). where Va = (B 2 /4iip) 1 / 2 is the Alfven speed, 
and B denotes the magnetic field strength. In a type II superconductor, B 2 is replaced by 
BH C \ in the formula for where H c \ ~ 10 15 G is the lower critical field, strengthening the 
coupling further. 

In a stratified star, the magnetized Ekman problem remains unsolved. Broadly, though, 
stratification (an d compressibility) inhibi t crust-core coupling, just as they do in an un- 
magnetized star ( lAbney fc Epstein! 119961 ). To estimate the stratified m agnetic couplin g 
time-scale, we follow the pioneering order-of-magnitude analysis in §6 of iMendelll ( 119981 ). 
which is arguably somewha t representative even though it applies to the idealized case 



of uniform magnetization. IMendelll ( 119981 ) showed that the t B formula in the previous 



paragraph applies only to the uppermost, buoyancy-limited boundary layer adjacent to the 
crust. One then imagines a stack of RQ/Va buoyancy-limited layers reaching from the 
crust to the core, which are activated sequentially and combine to give a revised time-scale, 
t B = 2(V A /Rty- 5/3 tt- 1 = 3.9 x 10 4 (p/po) 5/6 (fi/10 2 s- 1 ) 2 / 3 (5/10 12 G)- 5 / 3 s. Hence it is pos- 
sible to have t B > t u — and hence maintain core super-rotation on the viscous time-scale, as 
in §|5]and §3] — in objects with moderately magnetized interiors and/or fast spins, satisfying 
(fi/10 12 G) 5 / 3 (fi/10 2 s- 1 )- 2 / 3 < 5.5 x 10- 4 iUp/Po) 11/6 (T/10 8 K)~ 5 / 3 . Such objects may be 
found squarely within the main pulsar population in Figure [1] (e.g. B < 10 10 G) or among 
the recycled millisecond pulsars, pro vided that the internal fi e lds in the l atter are actually 
reduce d rather than just buried, cf. iKonar fc Bhattacharyal (119971 ) and iPayne fc Melatos 
( 120041 ) . The above conclusions are consistent with related work on boundary-layer damping 
of r-modes. In unmagnetized r-modes, it is found that compressibility approximately halves 
the damping rate, while stratification has almost no effect, because the radial fluid velocity 
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component in an r-mode is tiny ( IGlampedakis fc Anderssonll2006l ). unlike in the spin-down 
problem considered here, where t he radial and angular v elocity components are of the same 
order and stratification matters (lAbney &: Epstein! Il996l ). In magnetized r- modes, analyzed 
for the idealized, unstratified problem of a radial magnetic field with uniform or cosine- 
latitude magnitude, Lorentz forces shorten the damping time-scale, offset partially by the 
superfluid mutual friction, bec ause damping occurs in a boundary laye r abutting the crust, 
where buoyancy is negligible ( jMendelll l2001t iKinney fc Mendelll 120031 ) , unlike in core spin 
down. 

In a stratified star with a realistic internal magnetic geometry, things are more uncertain. 
Arguably, however, the overall impact of a nontrivial geometry is to lengthen the magnetic 
coupling time-scale further. One can see this qualitatively by considering the dispersion re- 
lations of hydromagnetic waves in a stratified slab, where the stratification axis n is inclined 
obliquely to both the background magnetic field B and the direction k of wave propaga- 
tion. A thorough treatment of this problem in plane-parallel geometry (without rotation , 
unfortunately) is given for various wave modes in Chapter 7 of iGoedbloed fc Poedtsl ( 12004 ) . 
The details fall outside the scope of this paper, but generically one discovers that the phase 
speed is typically lower than the buoyancy-limited Alfven speed in a uniform magnetic field 
and is much lower for certain orientations. For example, when n is perpendicular to k and 
Bo, and the wave propagates obliquely with respect to B , the phase speeds of the Alfven 
and slow modes (which are favored because they displace fluid horizontally, unlike the fast 
mode) tend to zero in the smal l-beta (pressure-dominated) regime relevant to neutron stars, 
as illustrated in Figure 7.10 of IGoedbloed &: Poedtsl ( 120041 ) . In the uniform-B calculations 
attempted so far in the neutron star literature, such orientations are avoided by construc- 
tion but t hey are inevitable in realistic geometries. For examp le, in a linked poloidal-toroidal 
structure ( jBraithwaite fc Nordlundl 120061 ; lLander et al.l 120121 ) , where field loops close inside 
the star, many internal surfaces exist where unfavorable orientations of n, k, and B are 
achieved, and buoyancy-modified hydromagnetic waves slow down dramatically, lengthening 
the crust-core magnetic couplin g time. In a tangled fi e ld, th is is even more likely to occur 
( jBraithwaite fe Nordlundl 120061 ). IGoedbloed fc Poedtsl ( 120041 ) emphasized in §7.3.3 of their 
book that these subtle effects emerge fully when oblique propagation (0 < |k ■ B | < 1) is 
considered, and that they are expected to be even more prevalent when n is inclined obliquely 
to Bp as well, the re l evant scenario in general in a neutron star [and outside the scope of 
Goedbloed fc Poedtsl (120041 )] . There is no substitute for a self-consistent, three-dimensional 
numerical simulation to quantify these difficult matters. The prospects are fair that suitable 
codes and computational platforms will become available shortly that are equal to the task. 



In a realistic magnetic geometry, another crucial complication arises: many field topolo- 
gies are incapable of enforcing crust-core corotation at all, no matter how fast the magnetic 
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coupling. Quoting from lEassonl (Il979al ). who demonstrated this counter-intuitive property 
analytically, "the field must have a topology such that the equation for the perturbation 
B' in the field induced by the slowing down has nonsingular, single- valued solutions". At 
first blush, this sounds like a side-issue of purely academic interest, but nothing could be 
further from the truth. There exist at least two important topologies (which are popular in 
the literature, incidentally) that violate the above condition and cannot enforce corotation: 
(1) axisymmetric, poloidal fields with loops closing inside the star; and (2) axisymmetric, 
toroidal fields. To see why, co nsider the mom entum equation in superfluid magnetohydro- 
dynamics, viz. equation (1) in lEassonl (Il979al ). Several of the forces are not conservative in 
general, e.g. the Lorentz force and the Fermi liquid interaction. @ The nonconservative forces 
cannot be balanced everywhere by pressure gradients and buoyancy, so B' tries to adjust 
until a balance is achieved. But often there is no way to do so. For example, if B is poloidal 
and axisymmetric, then B' is toroidal and the toroidal ((p-) component of the momentum 
equation reduces to (47r)~ 1 (B • W)(wB') = F(q,tu) in equilibrium, where w denotes the 
cylindrical radius, and the right-hand side is a complicated nonconservative function F(. . . ) 
of unperturbed ya riables q like the neutron and proton mass densities, viz. equation (6) 



m 



Eassonl (Il979al ). Integrating along a poloidal field line that closes inside the star, one 



reaches the unphysical conclusion that wB' is multi-valued. It is tempting to presume that 
the plasma can always find some way to adjust slightly and quasi-statically near a sensible 
topology to fix the pr oblem, " so that the relevant integrals along all closed loops give zero" 
(lEassonl Il979al ). But lEassonl (ll979af ) proved directly from the hydromagnetic equations of 
motion that the adjustment fails to achieve balance except in very special situations, because 
F(q, vd) is a function of variables that depend only on the unperturbed state of the star (den- 
sity, temperature, B , and the crust angular velocity, which do not respond as B' adjusts) as 
well as the angular velocity of the superfluid vortex array, which is governed by a superfluid 
equation of motion independent of B'. These important points are not appreciated widely 
in the literature and deserve careful study in future. For example, is the set of fields B that 
give unphysical solutions for B' "large" or "small"? lEassonl (Il979al ) expressed a hope that 
numerical simulations would answer this question, but the problem is formidable, and the 
relevant simulations have not yet been done. 

Magnetic flux is concentrated into quantized flux tub es in a type II supercond uctor, 
modifying the magnetohydrodynamic equations of motion (IGlampedakis et al.ll201ll ). The 



7 The entrainment force looks superficially like it is constructed out of gradients of the internal energy 
and proton effective mass, but the prefactors multiplying these gradients, like the proton number density 
and p ressure, are fun ctions of position even in equilibrium due to compositional stratification; see equation 
(2) in lEassonl |l979al) . 
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implications for core super-rotation remain highly uncertain. One outcome is to replace B 2 
by BH c i in the classic formula for the magnetic coupling time-scale, as described at the start 
of this section, but even this result has not yet been confirmed in realistic geometries and in 
the presence of stratification, where other length-scales and directional effects may enter the 
problem. Another outcome is that the flux tubes are expected to interact with the superfluid 
vortices in the neutron condensate, which are magnetized by entrainment. Opinions remain 
divided as t o whether the interact ion is strong, locking together the charged and neutral 
components ( Ruderman et al.l 19981) , or weak, as in the recent "snowplow" mo del of glitche s 



( IHaskell et al 



2012; 



Seveso et al. 



contains hyperons (jBabaev 



Andersson et al.l 120071 ; 



Link 



20121 ). if the core is a type I superco nductor (IJones 12006) 



2009), or is the seat of vortex instabilities ( IPeralta et al 



2006b 



20121 ). We do not attempt to adjudicate these difficult matters 



in this paper. Efforts are currently under way to calculate analytically the Ekman coupling 
time-scale in a type II superconductor including some of the above physics, albeit in an 
idealized plane-parallel geometry and a uniform magnetic field. ^ 

Whether strong or weak, magnetic coupling does not engage the uncondensed neutrons, 
which are viscous, participate in Ekman pumping, carry inertia, and therefore play a role 
in maintainin g core super-rotation [and in explaining radio timing data of pulsar glitch 



recoveries; see Ivan Eysden &; Melatosl ( 120101 )]. It is sometimes assumed that the uncondensed 
neutrons carry a tiny fraction of the star's inertia, because their density is proportional to 
exp(— A/hsT) in the ideal Bardeen-Cooper-Schrieffer t heory, and k^T is mu ch smaller than 
the superfluid gap A for most neutron pairing schemes ( lYakovlev et al.lll999f ). However, the 
ideal exponential gap formula tells only part of the story. Applied to superfluid helium, it 
predicts zero uncondensed atoms as T — > 0, whereas the experimentally measured fraction 
is ~ 14% due to van der Waals forces. It is likely that the nuclear color force is nonideal in 
the same way, as the quarks and gluons inside nucleons are polarized by random motions of 
nearby nucleons at high densities (jWeberl 120051 ). The effect has not been quantified yet from 
first principles, as it requires a hard calculation in quantum chromodynamics, but there are 
good indications that a van der Wa als equation of state is con sistent with measurements 
of the nucle ar liquid-gas trans ition (ITorrieri &; Mishustinl 120101 ) and also with scattering 
experiments (IKaiser et al.lll997l ). If the uncondensed fraction is substantial, it helps maintain 



residual core super-rotation via slow, unmagnetized, stratification-limited Ekman pumping, 
even when the other stellar components couple magnetically on a fast time-scale 



i pi 

.8 



8 K. Glampedakis, private communication. 

9 Recently an analytic theory of unmagnetized superfluid spin down has been developed in arbitrary 
geometry which self-consistently incorporates Ekman pumping, mutual friction, and the back reaction of 
the fluid on its container. The theory reproduces the results of classic laboratory experiments on liquid 
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5. Observational signatures 



How might one observe a crust-core shear of the magnitude in Figure HP The question 
is too multi-faceted to answer properly all at once, but we tender a few suggestions. Tak- 
ing i nspiration from the Ear t h, where seismic data carry the imprin t of the super-rotating 
core (jSong &: Richards! Il996l ; IZhang et al.l 120051 ; IWaszek et al.ll201ll ). one option is to look 
for rotational splitting of the quasi-periodic oscillat ions observed in the X-ray light curve in 
the immediate aftermath of a magnetar giant flare (jStrohmayer fc Watts! 120061 ) . At present, 
there is a lively debate as to whether magnetar oscillations arise from shear modes in the 
crust, magneto-inertial m odes in the core, a combination of the two, or something else en- 
tirely (jLander et al.ll2010l ). Theoretical treatments invariably start by perturbing a uniformly 
rotating background structure, and the rotational splitting is assumed to be small, as mag- 
netars are slow rotators. However, in objects where at least some of the natal rotation of 
the core is preserved as in Figure (U one might hope to see greater rotational splitting. 

Another possible signature is precession. A body composed of a rigid crust and liq- 
uid mantle precesses in a superposition of two modes (four if it also contains a solid in- 
ner core) excited by differen tial rotation, namely a Chandler wobble and free-core nutation 
(jKitiashvili fc Gusevll2008al ); think of an egg, whose yolk rotates differentially. The mode 
frequencies are related to the inertia of the components and t he strength of crust-co re cou- 
pling via magneto-inertial ( Baym et al. 1976 ) or elasto-inertial ( Bastrukov et ali ll999h 



waves. 



The model can be matched to the harmonics observed in PSR B1828— 11 (jStairs et al.ll2000 
Kitiashvili fc Gusev 2008a) jj] and to similar peri odic phenomena ranging from months to 



years in other objects (jKitiashvili &: Gusevll2008bl ) to infer the crust-core shear in principle. 
A long-period crust-mantle wobble may also be excited by a r otational glitch, when a portion 
of the crust-core shear is impulsively reset (jYuan et al.ll2010f ). 



If verified observationally, core super-rotation enriches neutron star astrophysics in two 
important ways. First, it implies ongoing, self- limiting, toroidal magnetic field generation. 
One might argue that magnetic stresses should react to erase the shear on the time-scale 



heliu m (jTsakadze fc TsakadzeHl98C| ) to 0.5% accuracy with zero free parameters |van Evsden fc Melatosl 



20111 ). yields excellent fits to Vela' s quasi-exponential glitch rec overies, and predicts without adjustment the 



"overshoot" observed in the Crab ( van Evsden fc Melatosl 2010 ). all with zero magnetic coupling. Of course, 
nonzero magnetic coupling is also consistent with the glitch data, provided that vortex pinning in the inner 
crust is invoked to retard the response; one is then led to postulate multiple pinning zones to match the 
multiple recovery time-scales observed in individual objects ( Alpar et al. 1984h . 

10 The oscillatory timing residuals and pulse profile changes in PS R B1828— 11 are now interpreted in 
terms of switching between two magnetospheric states (|Lvne et al.ll2010| ). 
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Ro _1 fi _1 , but the Sun provides a clear counter-example. Many resistive and other processes 
intervene to pinch off toroidal magnetic flux tubes before the back-reaction becomes too 
strong, e.g. the current-induced Tayler instability, which is especially active for high mag - 
netic Prandtl numbers (and is modified by the Hall effect in magnetars) (jRudiger et al.ll2009l ). 
or magnetic self-organization into a torus- like configuration, where Q i s constant along field 
lines, winding stops, and the magnetorotational instability is avoided (IDuez et al.ll2006l ). In 
numerical simulations of protoneutron stars, it is standard to find that the magnetic back 



reaction saturates at < 1% of the mechanical shear stress ( Cook et al.1 12003 



DuezetaL 



20061 ). Whether core super-rotation also drives a turbulent, self-sustaining, hydromagnetic 



dynamo in the face of strong stratification is less clear. Shear-driven turbulence in a strat- 
i fied spherical shell automatically supplies the non-zero helicity which a dynamo requires 
( IGellert et al.ll201lf ). and the magnetic Prandtl number easily exceeds the dynamo threshold 
~ 10 2 E 3 ^ in a neutron star. However, the interplay between S and a is not understood fully; 
recent work suggests that the dynamo depends weakly on S and generates a predo minantly 
dipolar field except for Ro ~ 1, but the jury is still out ( Simkanin fc Hejda 201 ll ). Two 
potentially observable manifestations are suggested by terrestrial analogies: geomagnetic 
jerks, wherein small changes in the length of the day are communicated to the magneti c field 
via a Hartmann current syst em connecting the Ekman/buoyancy layer to the interior ((Davis 



19941 ; iDavis &: Whaler! 119971 ). and longitud inal drift of equatorial geomagnetic spots, driven 
by zo nal flows over ~ 10 7 rotation periods ([Stanley &: Mohammadill2008l ; ISimkanin &: Hejda 



201 ll ). Spot or pole drift occurs over roughly th e same time-scale as the magnetospheric 



switching observed in objects like PSR B1828— 11 ( ILyne et al.ll2010l ). 

Secondly, core super-rotation opens up new channels of gravitational-wave emission 
( van Eysden fc Melatod l2008t Bennett et al.1 120101 ). The free energy in the sheared flow can 
excite nonaxisymmetric instabilities, e.g. the herringbone, Taylor-Gortler, and Stewartson 



11 iBraithwaitel fcood ) demonstrated simultaneously the twin phenomena of self-limiting toroidal field 
generation and a self-sustaining dynamo loop within a unified, self-consistent simulation of a differentially 
rotating, nonconvective star with and without stable stratification. Tayler and magnetorotational instabilities 
disrupt the toroidal field well before it grows to cquipartition, provided that the magnetic diffusion length- 
scale is shorter than the buoyancy length-scale (always true in a neutron star). Counter-intuitively, in 
the stratified system, a self-sustaining dynamo establishes itself more easily, the toroidal field builds up to 
saturation more slowly, and at saturation the toroidal field fluctuates less. In a statistically steady state, the 
magnetic back-reaction saturates at « 5% of the average magnetic pressure, much less than the mechanical 
stress driving the differential rotation. 



12 



Loperl (|197ll ) argued that magnetic and buoyancy forces ca n conspire to mimic unstratified Ekman 



pumping under special circumstances, clouding the picture further. Jiefu ( 19951) agreed but showed that the 
quasi-steady, sheared flow in the interior persists much longer than tg. 
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layer modes famil i ar from Newtonian and sup e rfluid spherical Couette fl ow (jNakabayashi et al 



20021 ; iHollerbachl I2OO3I: IPera 



counterparts (jSmirnov et al 



ta et al. 



2005 



2006a 



Peralta &: Melatosl 120091 ) and their stratified 



Munro et al. 2010 )^ l 13 l IH The theoretical machinery 



exists to treat such instabilities in a neutron star context, including stratificati on, com- 
press i bility, magnetic fields, and a dif ferentially rotating background in general (iMendell 
200ll ; iGlampedakis fc Anderssonl 120061 ). but it has only been applied to r- modes on a uni- 
formly rotating background at the time of writing. Core super-rotation also complicates 
gravitational-wave detection: the gravitational-wave phase may drift relative to the rota- 
tional phase inferred from radio pulse timing. Finally, as the Reynold s number E -1 is 



large, the interior flow is like ly to be turbulent, even when S is large ( iGreensteinl 11970 



Melatos fc Peraltal 120071 . |2010| ). Stratified turbulence collapses down to two dimensions (i.e 
spherical shells) only when the activity parameter Ro 3 5' _1 £' _1 drops belo w ~ 7, whi c h con- 



dition is avoided for Ro > 1 2 , i.e. throughout much of Figure [TJ [see llida et al.l (120091 ) 



and §2 in iMelatos fc Peraltal ( 120101 )] . Stochastic gravitational radiation from hydrodynamic 
turbulence imposes a fundamental noise floor on gravitational-wave observations of fossil 
rotators with Ro > 1 e ven with the current generation of long-baseline interferometers 



(IMelatos fc Peraltal 120101 1 . 



In summary, stratification limits Ekman pumping to a thin buoyancy layer in a neutron 
star, leaving a super-rotating core as the crust brakes electromagnetically. The crust-core 
shear grows, until dissipation (viscous or Eddington- Sweet) kicks in, then diminishes. Of the 
radio pulsars discovered to date, a handful of young objects (some Crab-like, some magnetar- 
like) are predicted to have unexpectedly high Rossby numbers Ro > 10 _1 , i.e. their cores 
are a fossil record of their angular velocities at birth. It is realistic to think that more 
candidate fossils will be discovered by future surveys. Moreover, there are 24 objects with 
Ro > 10~ 2 already known, some dominated by viscous diffusion, others by thermally driven 
circulation, which should be examined as a group for common features. If core super-rotation 
is confirmed, it will prompt a review of traditional ideas concerning how magnetic fields and 
gravitational waves are generated by neutron stars. 



This research was supported by a CSIRO-University of Melbourne Collaborative Grant. 



13 For example, iPeralta fc Melatos found Ro > AE 0A from numerical simulations for superfluid 

Stewartson layer instabilities, which is satisfied throughout much of Figure [TJ 

1 Nonax isymmetric flows exc ited in the core at birth should also persist up to the dissipation time-scale 
vam{t v ,t K ) ( Bennett et al. 2010h . 



Stratification is unimportant for r-mode damping, because r-mode circulation is nearly nonradial, unl ike 
spin down. Compressibility roughly doubles the r-mode damping time ( Glampedakis fe Andersson 20061 ). 
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